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Abstract 

The group velocity of light has been measured at eight different wavelengths between 
385 nm and 532 nm in the Mediterranean Sea at a depth of about 2.2 km with the 
ANTARES optical beacon systems. A parametrisation of the dependence of the 
refractive index on wavelength based on the salinity, pressure and temperature of 
the sea water at the ANTARES site is in good agreement with these measurements. 

Key words: ANTARES, neutrino telescope, optical beacon system, velocity of 
light, refractive index 
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1 Introduction 



The ANTARES neutrino telescope is located on the bottom of the Mediter- 
ranean Sea (42°50' N, 6°10' E) at a depth of 2475 m, roughly 40 km offshore 
from Toulon in France. The main objective of the experiment is the observa- 
tion of neutrinos of cosmic origin in the southern hemisphere sky. Sea water 
is used as the detection medium of the Cherenkov light induced by relativis- 
tic charged particles resulting from the interaction of neutrinos. The particle 
trajectory is reconstructed from the measured arrival times of the detected 
photons. The detector consists of 885 photomultiplier tubes (PMTs) mounted 
on twelve vertical lines with a length of about 450 m. The horizontal separation 
between hues is about 70 m. Further details can be found elsewhere [1,2,3]. 

Charged particles traveling through sea water produce the emission of Cherenkov 
light whenever the velocity of the particle exceeds that of light in water. The 
Cherenkov photons are emitted at a characteristic angle, 9c, with respect to 
the particle direction. This angle is related to the index of refraction of the 
medium as cos 9c = -j^- In this, /3 is the velocity of the particle relative to 
the speed of light in vacuum. The index of refraction, n^, corresponds to the 
ratio between the speed of light in vacuum and the phase velocity of light 
in water. The individual photons then travel through the water at the group 
velocity. Both the phase and the group velocity depend on the wavelength of 
the photons. This is usually referred to as chromatic dispersion. The group 
velocity is related to its phase velocity in the following way: 

where A is the wavelength of light. The index of refraction, Ug, corresponds to 
the ratio between the speed of light in vacuum and the group velocity of light 
in water. 

Since the PMTs cannot distinguish the photon wavelength, the variation of 
the photon emission angle and the group velocity due to chromatic dispersion 
cannot be accounted for on the individual photon level. Nevertheless, the av- 
erage effect of the wavelength dependencies are accounted for in the algorithm 
used to reconstruct the particle trajectory [4,5]. 

A measurement of the group velocity of light has been made using the optical 
beacon system of ANTARES. This system consists of a set of pulsed light 
sources (LEDs and lasers) which are distributed throughout the detector and 
illuminate the PMTs with short duration flashes of light. The refractive index 
is deduced from the recorded time of flight distributions of photons at different 
distances from the sources for eight different wavelengths between 385 nm and 
532 nm. 
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Fig. 1. Picture of a standard LED optical beacon (left), the modified LED optical 
beacon (middle) and a laser beacon (right). 

2 Experimental Setup 



The PMTs of ANTARES are sensitive to photons in the wavelength range 
between 300 nm and 600 nm. The maximum quantum efficiency is about 
22 % between 350 nm and 450 nm. The arrival time and integrated charge of 
the analogue pulse from the PMT are measured by the readout electronics [6]. 
The transit time spread of single photo-electrons of the PMT is around 3.5 ns 
(FWHM) [7]. 

The group velocity of light has been measured using the ANTARES optical 
beacon system. This system was primarily designed to perform time calibra- 
tion in situ [7,8]. There are two types of optical beacons, the LED optical bea- 
cons and the laser beacons. There are four LED optical beacons distributed 
along each detector line and two laser beacons at the bottom of two central 
lines. The in situ measurement of the temperature and salinity is provided by 
some Conductivity/Temperature/Depth sensor^. 

A standard LED optical beacon contains 36 individual LEDs distributed over 
six vertical faces forming a hexagonal cylinder housed in a pressure resistant 
glass enclosure (Figure 1, left). On each face, five LEDs point radially outwards 
and one upwards. All LEDs emit light at an average wavelength of 469 nm, 
except the two LEDs located on the lowest LED optical beacon of Line 12 
which emit light at an average wavelength of 400 nm. A modified LED optical 
beacon was installed in 2010 on Line 6. This LED optical beacon has three 
LEDs per face instead of six, all of them pointing upwards (Figure 1, centre). 
The three LEDs of each face emit light of the same colour. The average wave- 
length of the light from the six faces are 385, 400, 447, 458, 494 and 518 nm. 
The LEDs emit light with a maximum intensity of about 160 pJ and a pulse 
width of about 4 ns FWHM. The intensity of the emitted light can be varied 
changing the voltage feeding the LEDs. 

The laser beacon consists of a Nd-YAG solid state laser (Figure 1, right). It 
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Fig. 2. (a) Wavelength spectra of the hght sources used in this work as measured 
with a cahbrated spectrometer. The spectrometer measurements are made in air. 
The data points have been smoothed and the highest value of each spectrum is set 
equal to one. On the top of the Figure the peak wavelengths are indicated in units 
of nm. (b) Simulated light spectra at a distance of 120 m in sea water with the 
highest value renormalized to one. The differences between the spectra are due to 
the variation of the absorption length as a function of the wavelength. 

can emit pulses of light with an intensity of about 1 /iJ and pulse width of 
about 0.8 ns FWHM. The average wavelength is 532 nm. This light is spread 
out by an optical diffuser, so that the light can reach the surrounding lines. 
During calibration runs, the LEDs and lasers flash at a frequency of 330 Hz. 
Further details about the optical beacon system can be found elsewhere [7,8,9]. 

The wavelength spectra of the light sources used for this analysis were mea- 
sured using a calibrated spectrometei0(see Figure 2a). The typical width of 
each spectrum is around 10 nm except for the green LED (518 nm), which is 
larger, and the laser (532 nm), which is much smaller. 

Due to the wavelength dependence of the absorption of light in water, the spec- 
tra change as a function of the distance travelled by the light. The expected 
wavelength distributions as a function of distance have been estimated by 
Monte Carlo simulations using the dependence of absorption length on wave- 
length given by Smith and Baker [10]. In Figure 2b the spectra at a distance of 
120 m are shown for the different light sources. In particular, absorption has a 
large effect for wavelengths above 500 nm. Notice that the distributions have 
been renormalised to unity in each peak and therefore the relative effect of 
absorption between sources is not observed. This renormalisation is performed 
in order to show the change in the shapes of the distributions, which is what 
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Fig. 3. Example of a distribution of arrival times observed at a distance of 100 m 
from an LED. The inset shows a zoom around the signal region. The solid curve 
corresponds to a fit function as described in the text and given in Equation 2. 

influences the velocity measurement. The evolution of the spectra is taken into 
account in the final results (Section 4), in particular the uncertainty assigned 
to the wavelengths has been taken to be the root mean square (RMS) of the 
wavelength distribution given by the simulation. 



3 Data Acquisition and Analysis 



In order to measure the optical properties of the deep sea water, designated 
data taking runs were performed using the optical beacon system. During 
these runs, one single LED located in the lowest optical beacon of a line and 
pointing upward was flashed. Only the signals recorded by the PMTs along 
the same line are used in the analysis. As a result, the line movements due to 
the sea currents can safely be ignored. 

The runs used in this analysis were taken between May 2008 and April 2011. 
Each run contains typically more than 100,000 light flashes. Each flash is 
detected by a small PMT inside the optical beacon. The time of the flash and 
the arrival times of the photons on the PMTs were recorded within a time 
window from 1500 ns before to 1500 ns after the flash. The integrated charge 
of the analogue pulses of the PMTs were also recorded. Only runs were used 
when the average rate of background light was below 100 kHz. 

In Figure 3, the distribution of the arrival times of photons on a PMT located 
100 m above the LED optical beacon (A = 469 nm) is shown. The time, 
t = ns, corresponds to the time of the flash. A clear peak att = 470 ns can be 
seen which corresponds to the shortest propagation time of the light. The tail 
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Fig. 4. (a) Arrival time as a function of the distance between the LED (A = 469 nm) 
and the PMT for one single run. The solid line corresponds to a fit of a linear function 
to the data (see text), (b) Distribution of the measured refractive index for a total 
of 42 runs. 

with late photons can be attributed to hght scattering. The flat background 
arises from the optical background due to decays and bioluminescence. 

A convolution of a Gaussian and an exponential distribution on top of a fiat 
background is fitted to the data. The Gaussian distribution reproduces the 
transit time spread of the PMTs, the duration of the light flash and the effect 
of the chromatic dispersion in water. The exponential distribution takes into 
account the effect of the scattering of photons in water. The fit function can 
be formulated as: 

where t is the arrival time of the photons. The fit parameters are the optical 
background, B, the signal strength, S, the mean, /i, and width, a, of the 
Gaussian distribution and the exponential decay constant, r. In Equation 2, 
erfc(t) is the complementary error function distribution. An example of the 
fit is shown in the inset of Figure 3. The fit is determined in the range from 
200 ns before the most populated bin and 20 ns after. The arrival time of the 
light flash at each PMT is estimated by the fitted mean value of the Gaussian 
distribution. 

An example of the measured arrival times as a function of the distance between 
the optical beacon and the PMT is shown in Figure 4a. A linear function 
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Wavelength [nni] 


400 


470 


530 


Source of uncertciiiity 


Variation in [%] 


ivieuiioti 


0.27 


0.13 


0.04 


Absorption length 


0.11 


0.05 


0.07 


Scattering length 


0.22 


0.34 


0.24 


Total systematic uncertainty [%] 


±0.37 


±0.37 


±0.25 



Table 1 

The systematic uncertainty of the refractive index measurement is estimated for 
three different wavelengths by varying the absorption length and the scattering length 
(see text). 

has been fitted to the data to extract the group velocity of the hght. In the 
fitting procedure, the distance between the hght source and the PMT has been 
restricted to the range between 50 m and 250 m in order to hmit the effects 
of time slewing [7] and optical background on the measurement. 

A Monte Carlo simulation of the response of the detector to LED flashes 
has been made. The analysis method was performed to validate the analysis 
procedure and to study the systematic effects due to the assumed light absorp- 
tion and scattering parameters. Table 1 shows the different contributions to 
the systematic uncertainty in the measurement of the refractive index. These 
contributions have been determined as follows. For the wavelength values of 
400 nm, 470 nm and 532 nm, the difference of the refractive index obtained 
after the fitting procedure with respect to the input refractive index is shown 
in the first row of Table 1. This variation, which has been computed for the 
default values of absorption and scattering lengths, is termed "Method". In 
the second row, the variation in percentage of the fitted value is shown when 
the absorption length changes by ±50 %. Finally, in the third row the varia- 
tion in percentage of the fitted value is shown when the scattering length is 
varied between 20 m and 70 m [11]. As can be seen, the latter is the largest 
contribution to the systematic uncertainty. Adding in quadrature these values 
one obtains a systematic error that varies from 0.25 % to 0.37 % depending 
on the wavelength. 



4 Determination of the Refractive Index 

Between May 2008 and March 2010, a total of 42 runs were taken using an 
LED with an average wavelength of 469 nm. Three different LED intensities 
were used. For a high, middle or low intensity run the range of distances 
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Fig. 5. Index of refraction corresponding to the group velocity of light as a function 
of the wavelength. Also shown are results from measurements made in [11]. The 
grey band shows the systematic uncertainty. The two solid lines correspond to a 
parametrisation of the index of refraction evaluated at a pressure of 200 atm (lower 
line) and 240 atm (upper line) (see text). 

between the optical beacon and the PMT used in the fit were 50 - 250 m, 
40 - 220 m and 10 - 130 m, respectively. The measured refractive index values 
of these runs are shown in Figure 4b. In addition to these runs, 14 runs using 
an LED with an average wavelength of 400 nm and 13 runs using an LED with 
an average wavelength of 532 nm were taken. Between November 2010 and 
April 2011 eight runs with a modified optical beacon were taken, extending 
the measurements with six additional wavelengths. The index of refraction is 
estimated at each wavelength by the mean of the distribution. The measured 
refractive indices with the systematic uncertainty are shown in Figure 5 and 
tabulated in Table 2. As mentioned in Section 2, the uncertainties in the 
wavelengths have been taken to be the RMS of the corresponding distribution 
at the middle of the distance ranges. The variation of the RMS values in this 
range with respect to the middle is ±2 nm. 

The velocity of light in sea water at a given wavelength depends on the tem- 
perature, the salinity and the pressure of the water, because the density of 
sea water depends on these variables. A parametrisation of the light velocity 
proposed by Quan and Fry [12] is based on data from Austin and Halikas [13]. 
This parametrisation was modified to incorporate a correction for pressure [11]. 
During the data taking period, the temperature and salinity were measured 
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Table 2 

Summary of the refractive index results for the 2008 - 2011 data shown in Figure 5. 
For the refractive index first the statistical uncertainties are shown. They are com- 
puted as the RMS of the measured refractive index for each wavelength divided 
by the square root of the number of runs. In addition, there is a total systematic 
uncertainty of ib 0.005. 

in situ at a depth of 2250 m. At an ambient temperature of T = 12.9 ± 0.1 °C 
and salinity of S = 38.48 ± 0.01 %o, the refractive index corresponding to the 
phase velocity as a function of wavelength is expressed as: 



V — 200 

np{X,p) = 1.32292 + (1.32394 - 1.32292) x ^^^_^^^ + 

16.2561 4382 1.1455 x 10^ , , 
— + 

where A is the wavelength (in units of nm) and p is the pressure (in units of 
atm). Using Equation 1, the result of this parametrisation can be compared 
to the measurements (see Figure 5). 

From the known variations of temperature, salinity and pressure, the refractive 
index for a particular wavelength and at a given depth can be determined with 
an accuracy of better than 4 x 10~^. The parametrisation is in good agreement 
with the measurements. 

As mentioned in Section 1, the PMTs are unable to distinguish the wavelength 
of the incoming photons, so the effect of this chromatic dependence can only 
be taken into account on average. The spread of the arrival time residuals with 
respect to the expected arrival time of a 460 nm photon have been computed 
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by means of a standalone Monte Carlo simulation using the phase velocity 
for the emission angle and the group velocity (as given by the Equations 3 
and 1) for the arrival time. This simulation indicates that the spread of the 
time residual is 0.6 ns at 10 m, 1.6 ns at 40 m, 2.7 ns at 100 m and 3.6 ns at 
200 m. The time uncertainty introduced by this spread is unavoidable and is 
taken into account in the ANTARES official simulation program [14,15]. Even 
though the exact influence of the medium depends on the particular Cherenkov 
photon (wavelength, distance to the hit PMT) and therefore requires a full 
simulation, a rough estimate of the average effect can be obtained assuming 
that a majority of hits are between 40 m and 100 m from the track, which gives 
a value of ~2 ns for the uncertainty introduced by the transmission of light 
in sea water, including chromatic dispersion. This value is to be compared 
with ~1.3 ns coming from the PMTs transit time spread and to ~1 ns from 
time calibration. Monte Carlo simulation studies show also that an additional 
uncertainty of ~2 ns on the spread of the time residual degrades the neutrino 
pointing accuracy about 0.1 degrees. 



5 Summary 

Using pulsed light sources with wavelengths between 385 nm and 532 nm the 
group velocity of light in sea water at the ANTARES site has been measured 
as a function of wavelength. The emission spectra determined in the labora- 
tory for the different pulsed sources have been used as input to a Monte Carlo 
simulation in order to correct for the effect of absorption on the corresponding 
velocity measurement. Except for two sources these corrections are in general 
small. Likewise, a Monte Carlo simulation has been used to evaluate the sys- 
tematic uncertainties and to check, that the procedure to obtain the speed of 
light is robust and unbiased. The results obtained for the dependence of the 
group refractive index on wavelength are in agreement with the parametrisa- 
tion as a function of salinity, pressure and temperature of sea water at the 
ANTARES site. 
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